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Ionosphere scintillation can effect satellite telecommunication up to Ku-band. Nighttime 
scintillation can be attributed to large-scale inhomogeneity in the F-region of the ionosphere 
predominantly between heights of 200 and 600 km. Daytime scintillation has been 
attributed to sporadic E. It can be thought of as occurring in tluee belts: equatorial, high- 
latitude, and mid-latitude, in order of severity. Equatorial scintillation occurs between 
magnetic latitudes ± 25 degrees, peaking near ±10 degrees. It commonly starts abruptly 
near 2000 local time and dies out shortly after midnight. There is a strong solar cycle 
dependence and a seasonal preference for the equinoxes, particularly the vernal one. 
Equatorial scintillation occurs more frequently on magnetically quiet than on magnetically 
disturbed days in most longitudes. At the peak of the sunspot cycle scintillation depths as 
great as 20 dB have been observed at L-band. 


1. Introduction. . . , , , , , 

This paper is designed to respond to the question "What additional problerns does one have 
due to equatorial scintillation when operating satellite systems with equatorial earth stations 
at L-band". With this objective in mind we review the intensity, geomorphology and 
temporal characteristics of equatorial scintillations, discuss the indices which experimenters 
use to measure ionosphere scintillation, and look at currently available models. As the 
paper is aimed at non-ionosphericists, an effort is made to keep it as simple and 
applications-oriented as possible. 


1,1 uavivgi v/uhu. t 

Ionospheric scintillation was first recognized in the mid-forties (Flock 1987). Interestingly 
enough, the early observations of scintillation were of signals from radio stars - direct 
analogs of the earth-space transmissions of interest now. During the IGY there was special 
attention paid to the equatorial ionosphere and a host of strange effects were discovered 
which had the same temporal characteristics which we now associate with equatorial 
scintillation. One of these was named the "Far Eastern Anomaly" (Bateman, et al. 1959) 
and is illustrated in the beautiful range-time record shown in Figure 1, taken by Jim Watts 
on a low-latitude ionospheric scatter circuit. Note the sudden commencement of the 
broadened trace shortly after sundown which gradually decreases after midnight. 
Equatorial "bubbles" (regions of plasma deficiency) were observed at Jicamarca, Peru to 
follow the same temporal characteristics and were observed to rise to elevations of over a 
thousand km (Kelley et al. 1981, Kelley 1989). These plasma depletions are believed to be 
aligned along the Earth's magnetic field, have dimensions of 100-200 km east-west and 
several thousand km north-south, and to move eastward with velocities of 50-200 m/s. 
These depletion regions are then the seat of irregularities with scale sizes of meters to tens 
of kilometers (Johnson 1988). 

Persons interested in high-latitude scintillation are referred to an excellent paper by Kersley 
et al. (1988) on the subject. 
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1.2 Geographic and Temporal Characteristics. 

The oft-reproduced diagram shown in Figure 2 (Basu et al. 1988) contains a wealth of 
information. Solar cycle dependence is illustrated by the different depths of fading 
indicated for the solar maximum (up to 20 dB), and solar minimum (as low as 1 dB). The 
geographic extent is shown to be at least 20 degrees north and south of the equator at solar 
maximum, with peaks on either side of the equator. The diurnal behavior is in this Figure. 
Scintillations start shortly after 2000 local time and decay shortly after midnight. 

An ambiguity in Figure 2, as in similar earlier representations (e.g. CCIR 1986), is that the 
identity (geographic, geomagnetic or magnetic) of the equator is left vague. The reason for 
this is presumably that as worldwide measurements are not available it is difficult to 
generalize one's data to apply worldwide. 

By geomagnetic coordinates one refers to the Earth's dipole field, as represented by a 
magnet placed at the center of the Earth. The total intensity F of the Earth's dipole field in 
gauss is given by 

F = (0.3 R 3 /r 3 )[ 1 + 3 sin 2 @] ,/2 (1) 

where R is the Earth's radius (6371 km), r is the distance from the center of the Earth, in 
the same units as R, and 0 is geomagnetic latitude. The axis (slightly time variant) is given 
for January 1988 as passing through 79.1 degrees N geographic latitude and 289.1 degrees 
East geographic longitude (Davies 1989). Geomagnetic latitude for the Earth is given in 
Figure 3. The dip I in the geomagnetic coordinate system is simply related to the 
geomagnetic latitude L by 

tan L = 0.5 tan I (2) 

By magnetic or dip latitude L' (Smith 1957) is meant the latitude one would obtain from the 
analogous expression to (2) if the true ionospheric dip I' is substituted in (2), 

L' = tan 1 (0 5 tan I') (3) 

Magnetic latitude for the Earth is given in Figure 4. The normal procedure for obtaining 
magnetic latitude is to use the smoothed value of the surface dip and enter that in equation 
(3). As one goes up in altitude the direction of the Earth's field will gradually approach the 
dipole value given by expressions (1) and (2). How rapidly that value is approached is 
illustrated in Figure 5. For altitudes below 3000 km the true dip is seen to be better 
represented by magnetic latitude. 


2. Observations. 

The most applicable observations are the extensive ones made at C-band (primarily at 4 
GHz) at Intelsat stations and those at L-band on the 1.5 GHz signal from Marisat. There is 
also a limited number of recordings taken in August 1980 of LI and L2 GPS signals (Rino 
etal. 1981). 

2.1 Amplitude Scintillation. 

Shown in Figure 6 (Fang 1988) are recordings made at Stanley, the Hong Kong Intelsat 
station, on March 3-4, 1988 (at the peak of the solar cycle) of two satellites: the Pacific 
Ocean Regional (POR) at 174 degrees E and the Indian Ocean Regional at 63 degrees E. 
Hong Kong was the only receiver location. It could receive POR and IOR directly at 3950 
MHz and it could receive the signal uplinked to POR and IOR from Sentosa, Singapore and 
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Si Racha, Thailand at 6 GHz and down-linked at 4 GHz. If one assumes (as does the 
author) that the main effect is on the down-link (at the lower frequency) then one is faced 
with an eastward movement of the disturbance at the time when the sub-solar point is 
moving westward. This result is in agreement with spaced-antenna measurements of 
smaller scale inhomogeneity (Johnson 1988) but the rate of travel exceeds the 200 m/s 
found to be maximum in that experiment. 

2.2 Phase Scintillation. 

Figure 7, taken from Rino et al. (1981), is an example of data recorded during twelve 
consecutive evenings starting August 21, 1980 at Kwajalein Atoll. (Note that this was not 
a seasonal peak.) Measurements of amplitude and phase were made during scintillation 
periods on both the LI (1575.42 MHz) and L2 (1227.6 MHz) signal. During periods of 
moderate scintillation the phase deviations on LI exceeded 50 radians while those on L2 
exceeded 75 radians. In the example shown in Figure 7 the satellite (GPS 8) rose from the 
southeast around 2250 local time and achieved its maximum elevation of 70 to 75 degrees 
around 0330. 

To quote from page 256 of the article "During seven of these evening data runs, the 
scintillation was severe enough that the GPS 5010 receiver could not maintain lock on the 
signal during the early portion of the pass. When the satellite was above about 30 degrees, 
however, the scintillation was not severe enough to cause the receiver to break lock..." 

2.3 Frequency spectrum. 

By frequency spectrum is meant the prevalence of frequency components in the oscillations 
of the received signal under scintillation conditions. An example of this is shown in Figure 
8 where it can be seen that the most dominant frequency is around 10 seconds (O.l Hz). 
The prevalence of higher frequencies decays roughly as f 3 (CCIR 1990). 


3. Dependencies. 

The term "correlations" might be used here, but "dependency" seems to reflect a peak at a 
certain time of day or season somewhat better. 

3.1 Solar Cycle. 

There is a clear correlation of strong scintillation with sunspot cycle as reflected in Figure 9 
(CCIR 1990). However, as pointed out by Aarons (1981) this does not mean for all 
months or all times of day. Aarons relates the maximum scintillation regions to the 
Appleton anomaly, i.e. the fact that the maximum ionization in the F-region of the 
ionosphere is not located at the equator but in low-latitude peaks on either side of it as is 
illustrated in Figure 2. 

3.2 Season. 

The strongest scintillation at periods of high sunspot number for stations located near the 
two geographical maxima in Figure 2 are in equinoctial months with the vernal equinox 
preferred over the autumnal (Fang 1981, Basu et al. 1998). 

3.3 Diurnal. 

The diurnal variation of equatorial scintillation shows a strong preference for the hours of 
2000 to 0100 local time as is well illustrated by Figures 1 and 2. 

3.4. Magnetic. 

There are two aspects to the magnetic dependence. That concerning the dependence on 
geographic, dipole or dip latitude has been discussed under 1.2. The suggestion is that dip 
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or magnetic latitude is the better choice of the three. The other aspect is whether times of 
magnetic quiet or disturbance are preferred. As shown by Aarons (1981) and others there 
is a preference for quiet days for strong scintillation. 


4. Scintillation Indices. 

The quantification of scintillation indices can be attributed to Briggs and Parkin (1963). 
Their treatment was restricted to amplitude scintillation. They reasoned that, inasmuch as 
what is actually observed in recordings of satellite signals is a recorder deflection 
proportional to the amplitude of the received wave, measures related to the mean deviation 
and root-mean-square deviation were logical. These measures they denoted by Si and S 2 
defined as follows: 


Sj = k |r-r| 

R 


normalized mean deviation 


(4) 


s 2 = i {(R-R2)) 1/2 

R 


normalized RMS deviation 


(5) 


where R is the received amplitude, R is the mean value and R - R is the mean deviation of 
the sample record. The values are then normalized. Hence, Si is the normalized mean 
deviation, and S 2 is a normalized standard deviation. 


The authors reasoned further that, inasmuch as satellite scintillation recordings measure 
power, R 2 , rather than amplitude, R, two parallel indices, S 3 and S 4 , should be defined: 


s 3 = =L |r 2 - r 2 ! 

R 2 

(6) 

R 2 

(7) 


In their discussion Briggs and Parkin (1963) use an index S as a measure of scintillation 
depth defined as: 


s2 ?-(r? 

( 8 ) 



For a Rayleigh distribution, applying to the limiting case 
the relationships between the indices are: 

where phase deviations are large 

Si = 0.42 S 4 

(9) 

S 2 = 0.52 S 4 

( 10 ) 

S 3 = 0.73 S 4 

( 11 ) 

CO 

III 

CO 

-t*. 

( 12 ) 


206 


Whitney (1969) proposed a scintillation index, S.I., defined as: 

3 J — Pmax ~ Pmin (13) 

P max + Pmin 

where P m ax is the power amplitude of the third peak down from the maximum measured, 
and P m i n is the third peak up from the minimum excursion in a 15 minute period of 
recording. Later usage (Yeh et al. 1981) takes the highest peak and minimum excursion 
instead of the third up and third down. 


5. Models. 

Computer modeling of equatorial scintillation was initially based on two data sets (Secan et 
al. 1993). The first was data collected from the Wideband satellite in the 1970s at Ancon, 
Peru, and Kwajalein Atoll (VHF, UHF, and L-band). The second was data collected from 
the Marisat satellite from Ascension Island (L-band) and Huancayo, Peru (VHF), and from 
the Fleetsat satellite from Manila (VHF). The Wideband satellite data included both 
intensity and phase (WBMOD) and has been gradually enhanced over the years (Secan et 
al. 1987) to its latest version (Secan et al. 1993): "The WBMOD uses a collection of 
empirically derived models to describe the global distribution and behavior of naturally 
occurring ionospheric irregularities, and a power-law phase-screen propagation model to 
calculate estimates of the level of intensity and phase scintillation that these irregularities 
would impose on a user-defined system and geometry." Details of this code can be 
obtained from James A. Secan, Northwest Research Associates, Inc. P.O. Box 3027, 
Belleview, Washington 98009. 
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Figure 1. The "Far Eastern" or evening signal anomaly; after 
Bateman et. al., 1959) 
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Figure 2. Illustrative picture of scintillation occurrence based 
on observations at L-band (1.6 GHz); after Basu et al. 1988. 
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Figure 3. Geomagnetic dipole field coordinate grid 
superimposed on a Mercator projection of the world; after 
Johnson 1961. 




MAGNETIC LATITUDE DERIVED FROM THE MAGNETIC DIP 
Source: DTM 580 and Hydrographic Office Mop of Dip for 1945 
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Figure 4. Magnetic latitude derived from the observed dip, 
epoch 1955; after Smith 1957. 





Geographic Longitude =0° 

Source -Hydrographic Office Map for 1945 and CIW 580 



Figure 5. Comparison of the surface magnetic equator with the 
geomagnetic equator as a function of altitude above the earth; 
after Smith 1957. 
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Figure 6. Records of C-band ionospheric scintillation on March 
3-4, 1979 from Hong Kong; after Fang 1981. 
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Figure 7. Example of strong L-band scintillation and decay on = 
GPS signals; after Rino et al. 1981. 
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Relative power spectral density (dB) 



Fluctuation frequency (Hz) 

Figure 8. Power spectral density estimates for a geostationary 
satellite (Intelsat 4) at 4 GHz.; after CCIR 1990. 


The scintillation event was observed during the evenings of 28*29 April 1977 at Taipei Earth Station 

A: 30 minutes before event onset 

B at the beginning 

C: 1 hour after 

D: 2 hours after 

E: 3 hours after 

F: 4 hours after 
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Monthly sunspot numbers 

Figure 9. Dependence of 4 GHz equatorial scintillation 

on sunspot number; after Fang 1980, Fang and Pontes 1981 

The bars and squares are the ranges of variations over a year tor amerent earners. 

A: 1975-1976, Hong-Kong and Bahrein, 15 carriers 

B: 1974, Longovilo, 1 carrier 

C: 1976-1977, Taipei, 2 carriers 

D: 1970-1971, 12 stations, > 50 carriers 

E: 1977-1978, Hong-Kong, 12 carriers 

F: 1978-1979, Hong-Kong, 10 carriers 

G: 1979-1980, Hong-Kong, 6 carriers 
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